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Instanton Calculus of Lifshitz Tails 
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For noninteracting particles moving in a Gaussian random potential, there exists a disagreement 
in the literature on the asymptotic expression for the density of states in the tail of the band. We 
resolve this discrepancy. Further we illuminate the physical facet of instantons appearing in replica 
and supersymmetric derivations with another derivation employing a Lagrange multiplier field. 



Pure crystals have band structures with definite gaps 
in the energy spectrum. Impurities attach Lifshitz tails to 
each band [1] , allowing quasiparticles to be trapped deep 
in the band gaps [2] . Building on work by Halpcrin and 
Lax [3], Zittartz and Langcr [4] developed a systematic 
method for obtaining an asymptotic expression for the 
density of states deep in the tail of the band 
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for noninteracting quasiparticles moving in a random po- 
tential (see also [5]). A decade later, Cardy revisited 
the problem using the replica trick [6]. An instanton 
yielded the same exponential factor e~ B ( E \ but zero- 
mode counting showed that the way the prefactor A(E) 
scales with E is different from the one presented in [4]. 
To confirm our understanding of these methods, it is im- 
portant to reconcile the dispute [7] . 

It turns out that the source of the disagreement lies in 
a minor algebraic mistake. Correcting Eq.(5.24) of [4] to 
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the apparent discrepancy between [4] and [6] is resolved. 
The eventual agreement adds confidence to the use of the 
replica trick in the nonperturbative regime. 

In the remainder of this note, we survey various meth- 
ods for obtaining the asymptotic forms of Lifshitz tails, 
illuminating connections among them and tying up loose 
ends. We first present a derivation utilizing the method 
of Lagrange multipliers, building upon classic work by 
Halperin, Lax, Zittartz, and Langer [8]. It adheres a 
physical interpretation to the instantons appearing in su- 
persymmetric and replica derivations, which we review 
next. Common to all three derivations is the appear- 
ance of translations! zero modes, but each method treats 
remaining modes differently. In particular the supersym- 
metric method replaces bosonic zero modes appear- 
ing in the replica limit by the combination of one bosonic 
zero mode and two fermionic zero modes. 

Focusing on the vicinity of a band edge, we consider 
a system of noninteracting quasiparticles in d spatial di- 
mensions, governed by the Schrodinger equation 

-|V<(x) + K(x)^(x) = i^(x) (3) 



where V(x) is a Gaussian random potential which is as- 
sumed to self-average on macroscopic scales. In the ther- 
modynamic limit where the volume of the sample Vd ap- 
proaches infinity, we can then legitimately estimate the 
disorder-averaged density of states as 
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Here, A/j = / [ W] e iyf dxV2 W is the normalization 
constant and 7 characterizes the strength of the disorder. 
We are interested in the asymptotic behavior of [p(£')] d . 
in the limit of large negative E. 

To this end, we seek a localizing potential which mini- 
mizes the cost j- J dxV 2 (x) while still holding an eigen- 
function with the negative eigenenergy E. Note that for 
d < 4 a square-integrable potential always has a unique 
normalizable ground state [10] and this is the eigenfunc- 
tion we seek henceforth as, otherwise, we can lower the 
cost of the potential via simple rescaling. Through the 
introduction of a Lagrange multiplier field A (x) and a 
Lagrange multiplier /io, the problem becomes equivalent 
to the minimization of the cost action 



I[V (x),V(x),A(x),H 
= + i/M 2 (x) 

~ J dxA(x)/i?+|V-F(x) 



(6) 



-Mo 



dxtp 2 (x) - 1 



Extremizing it yields 

7(x) = -A(x)V(x) ; 
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V 2 4' (x) + V (x) ijj (x) = Eip (x) , 
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/ dx-0 2 (x) = 1, fio = 0, and 

r> 2 



— V 2 A (x) + V (x) A (x) = EX (x) 
2 m 
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The last equality, combined with the uniqueness of the 
ground state, dictates that A (x) = AoV' (x) with a con- 
stant Ao- Our problem then morphs into the instanton 
problem with a single real scalar field, which we know 
it to have spherically symmetric solutions minimizing 
the action among all the nontrivial stationary points for 
d < 4 [11]. Thus we have the cost minimizing solutions 
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where /(f) satisfies 
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ar z r dr 
= 0, monotonically decreasing toward as 



f — > oo. The normalization condition on ipo fixes Ao 
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Evaluating the cost action for these solutions yields the 
leading exponential factor e~~ with the dimensionless 
coupling 



and the dimensionless number 

r°° 
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To evaluate the sublcading prcfactor, let us first ex- 
pand around each saddle as 



V-V = J2& v i (16) 

1=0 

where v^s are a set of orthonormal functions. We choose 
v = A E f 
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with A E = (2a d )-*(-E)* (|£) 4 so that all the other 
modes will not change the ground state energy to first 
order in Integration over £o, hitting the energy delta 
function in the density of states, leaves us with the factor 



m — dE 
d£, 




(18) 



where we also took care of the factor coming from A/j 
Next for i = 1, d we choose 



Vi = A T diV 



where A T = c T {-E)i-i (|£) 4 2 
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with 
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They are d translational zero modes and integration over 
these modes should be traded for integration over the 
collective coordinates xo, sweeping along the saddle sub- 
manifold in the field space. The Jacobian involved in 
this coordinate transformation is A^} for each mode as 
can be seen by comparing changes in the field induced by 
(i5xo) i and by After dividing by the volume Vd and 
again taking A/j into account, we receive 
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from these modes. 

Finally integration over all the other modes gives a 
constant of order 1, worked out in the Appendix. 

All in all we find 



[P(%,^x(-£)i 
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{g(E)} 2 xe f< E > 
(22) 
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for d < 4 in the regime E <C —73=3 (2™.) i - d . An ex- 
pression for the dimensionless overall constant c is given 
in the Appendix [cf. Eq.(47)]. 

Let us now turn to the supersymmetric derivation [15— 
17]. First we express the density of state as 



p v (E) = — lim Im 

7T 5^+0 



i / dxG^(x,x;E + i5) 
Vd J 
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The retarded one-particle Green function Gp(x,x';i? 
i5) can be represented as 



/ [V<f>} 0(x)</-(x')e lSl 
$[D(f\e iS v\ 
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The supersymmetric method proceeds by rewriting the 
expression (24) as 



0(x)-0(x )e 
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with 



and 



(27) 

where we doubled the bosonic field <j) to (f> = (<pi, 4>2 ) and 
introduced fcrmionic fields xi and Xi- Now that there is 
no denominator containing the random potential, we can 
perform functional integration over V and obtain 
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, (29) 
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S2i[xi,X2,<fi] = ^ dx X2 i--^^ 2 - E - -<j) J 



and 
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(30) 
(31) 



We have defined <fi = e + ~<p and Xa = e + ~Xa for E < 
and the expression (28) should be viewed with appropri- 
ate analytic continuation in mind [18, 19]. 

To evaluate [p{E)] d a for large negative E, we use the 
method of steepest descent, extremizing Sh- The trivial 

saddle <f> = gives no contribution to [p(£/)] d due to the 
absence of negative modes. Among nontrivial saddles, we 
assume [20] that the saddles 




x - x 
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Here, e_j_ is a unit vector perpendicular to e, tpf s are a 
set of orthonormal functions satisfying 



.^-E-^f\^={-E)4^\ (36) 



and (pf 1 - = ipf = s are another set of orthonormal 



functions satisfying 



with dimensionless numbers cl s and c^'s. We deal first 

B I 

with £ ; fluctuations and then with the rest. 

B I 

Analyzing (p, modes, we find that the lowest mode 
has a negative eigenvalue cjj < 0, giving rise to a factor 
of i(—E)~i and allowing the saddles to contribute to the 
density of states. Next come d translational zero modes. 
Trading integration over these modes for integration over 
Xq and dividing by the volume, we receive the Jacobian 



rfx V(f> c 
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(38) 

Finally all the other modes have positive eigenvalues, 
each of which gives a factor of (— E)~^. 

Analyzing the other set of fluctuations, except the 
lowest modes, all the modes have positive eigenvalues, 
each of which gives a factor of (— _E) _ 2 + 2+3 = (— £ , )+3. 
The lowest modes are the zero modes arising from 0(2)- 
rotational symmetry, proportional to <f> cl (x) . The 

bosonic zero mode, upon trading integration over 
for integration over e, yields the Jacobian 



minimize the action, where e is a constant unit vector 
and /(f) was defined around Eq.(12). Evaluating the ac- 
tion for these solutions gives the same leading exponential 
factor e~~^ as before. 

In regards to the subleading pref actor, one contribu- 
tion comes from 



dx</> cl (x) • </> cl (x) ~ g 1 (-E)- 



(33) 



in front. (We will not keep track of the overall dimen- 
sionless constant in this derivation.) To evaluate the re- 
maining contributions, we expand around each saddle as 



EtB|| B|| -. V^B_L,_B_L 



(34) 
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dx.<p cl 
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We also need to saturate fermionic zero modes by ex- 
panding the action to the quartic order in fluctua- 
tions: if we kept only quadratic terms in the expan- 
sion of the action, integration over (£o) a ' s would give 
zero. Thus we must bring down either a factor of 



lfdxx2 N> -</>oiJ Xi or jfdx(x2Xi) ■ After appro- 
priate Gaussian integrations, we obtain a factor of 



(-EY 



(40) 



Putting them all together, we recover the same result 
(22) as before. 



4 



Finally let us turn to the replica derivation [6]. The 
replica trick proceeds by rewriting the expression (24) as 



so as to keep the ground state energy intact to this order. 
The resulting disorder cost is 



—i 



(x) .^( x ') e ^f=iSv[ 



(41) 



where we introduced (N r — 1) replicas, promoting <\> to 
4> = (4>i, 4>2i ■■■■> 0JV r )j and took the dicey limit in which 
N T —> to eliminate the denominator. After integrat- 
ing over V and making analytic continuation, we find 
instantons of the same form (32), but zero-mode analysis 
is slightly different from the one in the supcrsymmetric 
derivation. Besides d translational zero modes, there are 
lini7v r _ i .o(-A r r — 1) = — 1 bosonic zero modes coming from 
0(-/V r )-rotational symmetry. The latter is replaced by the 
combination of one bosonic zero mode and two fermionic 
zero modes in the supersymmetric derivation. 

We note that the instantons (32) appearing in replica 
and supcrsymmetric derivations and the localized wave- 
functions (11) have exactly the same shape. Thus we 
interpret the instantons as most likely forms of localized 
wavefunctions or square roots of localizing potentials [cf. 
Eq.(lO)], dilutely distributed for large negative E. It may 
be more appropriate to call all these solutions "localons." 

Through the use of the Lagrange multiplier field, we 
simplified the method developed by Halperin, Lax, Zit- 
tartz, and Langer. Instantons appearing in replica and 
supersymmetric derivations are now interpreted as local- 
ized wavefunctions or as square roots of localizing poten- 
tials. All three derivations are shown to give the same 
asymptotic expression for the Lifshitz tail in the case of 
the Gaussian random potential, including the subleading 
prcfactor. 

The author thanks Steven A. Kivelson for pointing out 
the discrepancy between [4] and [6] and for emphasizing 
the conceptual importance of the would-be conflict. He 
thanks Allan W. Adams, John L. Cardy, John A. Mc- 
Greevy, and Stephen H. Shenker for comments, discus- 
sions, and encouragement. He is supported by a JSPS 
Postdoctoral Fellowship for Research Abroad. 



APPENDIX 

For I > d + 1 , we have the ground state energy shift 
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to second order in £/'s. We compensate it by setting 
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Imitating [4], we proceed by choosing v\ = fui with 

h 2 



(44) 
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V 2 - E + (1 + Cl )V 



ui = 0. 



(45) 



(Corresponding to vo and Vi, we have uo oc / and m oc 
dif with Co = and c, = 2, respectively.) With this trick 
we evaluate the cost to be 
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Performing Gaussian integrals over £;'s for I > d + 1, 
taking A/j into account, and combining with the contri- 
butions from the other modes, we obtain 



. 27T2( 
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l=d+l 



(47) 



Another expression for c is given in [9] . 

For d = 1, we get = |, c\ = 4, ct = ^p-, and 
c; = the last of which can be obtained through 

the use of Gegenbauer polynomials of order | [4]. Thus 
c = — , conforming with the exact result obtained by 
Halperin [22]. 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 



[8] 

[9] 
[10] 

[11] 



I. M. Lifshitz, Adv. Phys. 13, 483 (1964). 

P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

B. I. Halperin and M. Lax, Phys. Rev. 148, 722 (1966). 

J. Zittartz and J. S. Langer, Phys. Rev. 148, 741 (1966). 

B. I. Halperin and M. Lax, Phys. Rev. 153, 802 (1967). 

J. L. Cardy, J. Phys. C: Solid State Phys. 11, L321 

(1978) . 

The author thanks Steven A. Kivelson for pointing out 

the conflict during his lecture course on disorder and for 

stressing the importance of the problem. The resolution 

of the conflict owes a great deal to him. 

We perform the restricted functional integral in a way 

different from the one presented in previous work [9]. 

There they introduce a Lagrange multiplier while here 

we introduce a Lagrange multiplier field. 

A. Houghton and L. Schafer, J. Phys. A A 12, 1309 

(1979) . 

E. H. Lieb and M. Loss, Analysis, 2nd ed. Graduate Stud- 
ies in Mathematics, Vol. 14, 2nd ed. (American Mathe- 
matical Society, Providence, RI 2001). 
The case for d — 3 is worked out in [12] and the case for 
d = 2 is worked out in [13, 14]. For d = 1, /(f) 
gives a desired solution [cf. Eq.(12)]. 



V2 
cosh(f) 



5 



[12] S. R. Coleman, V. Glaser and A. Martin, Commun. 

Math. Phys. 58, 211 (1978). 
[13] H. Berestycki, T. Gallouet and O. Kavian, C. R. Acad. 

Sci. Paris 297, 307 (1983). 
[14] H. Brezis and E. H. Lieb, Commun. Math. Phys. 96, 97 

(1984). 

[15] J. M. Luttinger and R. Waxier, Ann. Phys. 175, 319 
(1987). 

[16] Th. M. Nieuwenhuizen, Physica A 167, 43 (1990). 

[17] For a review of the supersymmetric method, see K. Efe- 

tov, Supersymmetry In Disorder And Chaos (Cambridge 

University Press, Cambridge, 1997). 
[18] J. S. Langer, Ann. Phys. 41, 108 (1967) [Ann. Phys. 281, 

941 (2000)]. 



[19] For a review, see S. Coleman, Aspects of Symmetry (Cam- 
bridge University Press, Cambridge, 1985). 

[20] For d = 3, it has been argued that the saddles of the 
form (32) have the least action among all the nontriv- 
ial stationary points [21], extending the theorem in [12]. 
For d < 2, the author is unaware of analogous extension 
to multiple fields needed in replica and supersymmetric 
derivations, but the resultant agreement with the other 
derivaton strongly indicates that they indeed minimize 
the action among all the nontrivial saddles. 

[21] J. M. Luttinger and R. Tao, Ann. Phys. 145, 185 (1983). 

[22] B. I Halperin, Phys. Rev. 139, A104 (1965). 



